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We investigate the propagation and scattering of liighly nonlinear waves in disordered granular 
chains composed of diatomic (two-mass) units of spheres that interact via Hertzian contact. Using 
ideas from statistical mechanics, we consider each diatomic unit to be a "spin" , so that a granular 
chain can be viewed as a spin chain composed of units that are each oriented in one of two possible 
ways. Experiments and numerical simulations both reveal the existence of two different mechanisms 
of wave propagation: In low-disorder chains, we observe the propagation of a solitary pulse with 
exponentially decaying amplitude. Beyond a critical level of disorder, the wave amplitude instead 
decays as a power law, and the wave transmission becomes insensitive to the level of disorder. We 
characterize the spatio-temporal structure of the wave in both propagation regimes and propose a 
simple theoretical interpretation for such a transition. Our investigation suggests that an elastic 
spin chain can be used as a model system to investigate the role of heterogeneities in the propagation 
of highly nonlinear waves. 

PACS numbers: 05.45.Yv, 43.25.-|-y, 45.70.-n, 75.10.-b,46.40.Cd 



I. INTRODUCTION 

A recent focal point of studies on nonlinear oscilla- 
tors is one-dimensional (ID) granular crystals, which con- 
sist of closely packed chains of elastically colliding par- 
ticles. Because of the contact interaction between parti- 
cles, granular crystals are characterized by a highly non- 
linear dynamic response. This has inspired numerous 
studies of the interplay between nonlinearity and dis- 
creteness [ll, Q ■ Granular crystals can be created from 
numerous material types and sizes, making their prop- 
erties extremely tunable Because they provide an 
experimental setting that is both tractable and flexible, 
granular crystals are an excellent testbed for investigat- 
ing the effects of structural and material heterogeneities 
on nonlinear wave dynamics. Recent studies have exam- 
ined the role of defects [a], interfaces between two dif- 
ferent types of particles 0,13: decorated and/or tapered 



chains [8|, |9| , chains of diatomic and triatomic units 10 - 
12 1, and quasiperiodic and random configurations [14- 
16|. The tunability of granular crystals is valuable not 
only for basic studies of the underlying physics but also 
in potential engineering applications, including shock and 
energy absorbing layers [a. Ha . H. sound focusing devices 
and delay lines 13|, actuators [isll. vibration absorption 
layers [ll| , and sound scramblers [3, [3 • 

An important context for the present paper is the con- 
sideration of disordered settings in nonlinear oscillator 
chains. Some of the most prominent recent investiga- 
tions in this area have generalized to weakly nonlinear 
settings [l^l the ideas of P. W. Anderson, who showed 
theoretically that the diffusion of waves is curtailed in 
linear media that contain sufficient randomness induced 
by defects or impurities [2l|, [2^ . Research on the An- 



derson transition is extremely exciting and is germane to 
the investigation of disordered systems more generally, 
but our goal in this study departs somewhat from this 
thread, as we seek to investigate order-disorder transi- 
tions in strongly nonlinear media such as granular crys- 
tals. This is a key challenge in the study of nonlinear 
chains [23| . 

In this paper, we investigate ID granular crystals 
that consist of chains of elastic spheres that interact via 
Hertzian contacts. The chains are composed of diatomic 
units: two-mass cells (henceforth called "dimers") that 
each consist of one particle of one type ("type 1") and 
one particle with different properties ("type 2"). Each 
dimer can be arranged in one of two orientations: type 
1 followed by type 2 or vice versa. To quantify the na- 
ture of the disorder and define the chain configuration, 
we borrow an idea from statistical physics and treat each 
cell as a "spin" (where, for example, spin "up" is charac- 
terized by a dimer composed of a type 1 particle followed 
by a type 2 particle, and spin "down" is characterized 
by a dimer composed of a type 2 particle followed by 
a type 1 particle). This allows a novel perspective on 
investigations of granular crystals, as one can examine 
the wave propagation dynamics as a function of an order 
parameter, which we define using the ratio of the num- 
ber of spins oriented in the same direction to the total 
number of spins. Borrowing another technique of statis- 
tical physics (which, in particular, is used in the context 
of disordered systems), we take an ensemble average of 
quantities that characterize the wave propagation (am- 
plitude, wave velocity, etc.) over chain configurations 
with the same level of disorder. Accordingly, as disorder 
increases, we are able to show that the wave transmission 
decreases rapidly, saturating to a value that is indepen- 



2 



dent of the chain heterogeneity. 

To characterize the order-disorder transition, we an- 
alyze local properties of the wave (namely, its ampli- 
tude and spatial structure) during its propagation. This 
reveals two propagation mechanisms: We find solitary 
waves below a threshold level of disorder and delocal- 
ized waves above it. We use the term "wave delocaliza- 
tion" to mean that the spatial structure of the solitary 
waves — which extend over only a few beads and are ob- 
served in the low-disorder regime — disappears with in- 
creasing disorder, giving rise to a significantly more ex- 
tended wave structure that is present over a long portion 
of the chain. We also investigate the effects of chain 
length and material properties and propose a simple the- 
oretical description — based on a scattering mechanism 
across individual defects in the chain — that accounts for 
most of the wave properties. 

The remainder of this paper is organized as follows. 
We detail our definition of elastic spin chains in Sectionllll 
and discuss the setup for our experiments and numerical 
computations, respectively, in Sections IIIII and IIVI We 
then discuss in more detail the order-disorder transition 
and the dynamics in the low-disorder and high-disorder 
regimes in Sections fVl and IVTl We discuss our results in 
Section IVlIl and conclude in Section IVIIIl 



II. ELASTIC SPINS AND DISORDER 

The concept of spin has been used successfully to de- 
scribe myriad physical phenomena, such as magnetiza- 
tion and glassiness [23|. In this paper, we apply this 
idea to analyze order versus disorder in granular crys- 
tals. Fully ordered diatomic chains allow the propagation 
of solitary waves [Toj . However, an orientation reversal 
of even one cell causes a defect in the system, leading 
to complicated dynamics such as partial wave reflections 
and radiation shedding ■ Increasing the hetero- 

geneity of the chain further via random arrangements of 
cells alters the dynamics of wave propagation drastically 
and necessitates a different approach. To understand 
such heterogeneous granular chains, we characterize the 
level of disorder in the chain by using the diagnostic 

jj ^ y I^UP ~ ^downl (^^^ 
^up + ^down 

where iV^p is the number of dimcrs (spins) with one 
orientation (called "up") and -^Jown number of 

dimers with the opposite spin (called "down"). The total 
number of spins, = iVup + ^down' equal to one half 
of the number of particles in the chain. By convention, 
we say that a dimer has spin "up" if the heavier particle 
is nearer the direction of the original incident wave. The 
definition ([!]) guarantees that maximum order — the case 
studied in Ref. [T^], in which all dimcrs have the same 
orientation — occurs when D = Q. Minimum order (and 
hence maximum disorder) occurs when A'up = -^down' 



which yields D = 1 [s^l- The quantity D thereby mea- 
sures the amount of heterogeneity in the granular chain. 
It is worth remarking that D = 2Ndefect/N, where iVdefect 
is the number of "defects" (situations in which consecu- 
tive spins have opposite orientation) in the chain. 

One obtains a given value of D for many different possi- 
ble combinations of dimcrs, so D indicates only the pres- 
ence of "defective" spins and not their location. To ac- 
count for this, we average over multiple disordered con- 
figurations with equal values of D. The effect of the 
level of disorder D on the transmission of nonlinear waves 
through the chain is the central point of this study. As 
D increases, we expect the amplitude of the transmitted 
force to decrease. However, our study reveals (see the 
discussion below) that this is true only for sufficiently 
small values of the parameter D, for which the solitary- 
like wave structure is preserved. In contrast, we observe 
waves with delocalizcd structure at larger values of D, 
and the wave propagation is no longer affected by the 
level of heterogeneity in the chain. The issues that we 
address below are (1) when this transition occurs and 
(2) how we can characterize and interpret both of these 
regimes. 



III. EXPERIMENTAL SETUP 

To experimentally investigate the order-disorder tran- 
sition in granular chains, we construct chains of dimer 
cells composed of units of spherical beads of two different 
materials and/or sizes. For the experimental results in 
our main discussion, we use 4.76 mm radius stainless steel 
(non-magnetic, 316 type) and aluminum (Al; 2017- T4 
type) spheres. The steel and Al beads, respectively, have 
masses 3.63 g and 1.26 g, elastic moduli 193 GPa and 
72.4 GPa, and Poisson ratios 0.30 and 0.33 [Ulill. We 
examine chains with various sequences of up and down 
cells and assemble each chain in a horizontal setup [see 
Fig. [Ha)] that is composed of two steel bars clamped on a 
sine plate. We ensure contact (and hence ID impact and 
wave propagation) between the particles by tilting the 
guide slightly (3.5 degrees). The chain is 31 cells long. 

We generate the incident solitary waves by impacting 
the chain with a steel particle "striker" (of the same size 
as the steel particles in the chain), which is launched 
down a ramp. Based on the drop height of the striker, we 
calculate its impact velocity to be about 0.5 m/s for these 
experiments. To visualize the traveling waves directly, we 
place piezo sensors {RC ^ 10'^ /iS, Piezo Systems Inc.) 
inside small (2.38 mm radius) stainless steel beads (non- 
magnetic, 316 type) [l^jIH,!!^] that we use in place of the 
light Al particle in the 20th and 25th cells. The sensors 
are kept in place by thin support plates with circular 
holes whose diameter is slightly larger than that of the 
beads. 

Using this setup, we consider values of the parameter 
D, which characterizes the level of disorder in the first 
20 pairs of the chain (see the discussion above), rang- 
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FIG. 1: [Color online] Schematic of the experimental setup. 



ing from to 1 in increments of 0.1. For each value, 
we consider 3 different dimer cell arrangements, and we 
average the results from 3 striker drops for each such 
arrangement. We also performed similar experiments 
with other types of diatomic chains — including one com- 
posed of large stainless steel (radius 4.76 mm) and small 
(2.38 mm) stainless steel spheres and another composed 
of small stainless steel and small PTFE (also 2.38 mm) 
spheres. We obtained consistent results with all of the 
considered configurations. In our discussion below, we 
present the steel-aluminum example. 



IV. NUMERICAL SETUP 

We model a chain of 2A'' spherical beads as a ID lattice 
with (conservative) Hertzian interactions between parti- 
cles [!|,[2: 



A 



ij + l r3/2 



1/2 



'j+i 



1 - i^f 



(2) 



where j £ {1, • • • , 2A^}, yj is the coordinate of the cen- 
ter of the jth particle measured from its equilibrium po- 
sition, 6j = max{j/j„i — yj,0} for j S {2, 3, . . . , 2A^}, 
5i = 0, 52N+1 = niax{y2Ar, 0}, Ej is the elastic modu- 
lus of the jth particle, Vj is its Poisson ratio, rrij is its 
mass, and Rj is its radius. The particle j = repre- 
sents the striker, and the (2iV + l)th particle represents 
the wall. The initial velocity of the striker is determined 
from experiment, and all other particles start at rest in 
their equilibrium positions. 



V. COMPARISON BETWEEN EXPERIMENTS 
AND NUMERICS 

We compare numerical simulations of Eq. ^ with ex- 
periments using the same type and quantity of particles 
composing the chain, the same "spin configuration" , and 
the same excitation input. Once the striker hits the first 
bead, a compression wave propagates through the chain. 



We measure the force of the propagating wave using a 
sensor that we placed in the 20th cell. In Fig. [DJa), we 
show plots of force versus time for both low and high 
values of D for the experiments and the numerical com- 
putations. Panel (b) of Fig. [2] shows the peak amplitude 
of the transmitted wave as a function of D for both ex- 
periments and numerical simulations. In both plots, we 
show the normalized force (which is given by the mea- 
sured force divided by the peak force that we obtained in 
the perfectly ordered case). We consider three arrange- 
ments corresponding to each D value, and we acquire 
three independent measurements for each arrangement 
to verify experimental reproducibility. In Fig. [2Ka), we 
show the force-time history averaged over all three dif- 
ferent dimer cell arrangements at each level of disorder. 
In Fig. HKb), we show the peak normalized transmitted 
force. Each color/shape corresponds to a given spin con- 
figuration at the stated value of D, and each marker rep- 
resents a single experimental run. The solid curve gives 
the median transmitted force value for each level of dis- 
order. The squares give the results of the numerical runs 
that correspond to the median experimental configura- 
tions. 

We observe good qualitative agreement between nu- 
merics and experiments. Note, in particular, the de- 
crease of the peak force as D becomes larger. Because the 
numerics and experiments agree qualitatively — the small 
quantitative difference arises from experimental dissipa- 
tion that is not modeled in Eq. ([2]) — we focus the remain- 
der of our discussion on the former. 



VI. ORDER-DISORDER TRANSITION: FROM 
SOLITARY TO DELOCALIZED WAVES 

A. Wave Transmission Model 

Experiments and simulations (with as many as A'^ > 
500 cells in the latter) both suggest the existence of two 
different regimes: At low disorder — i.e., for D smaller 
than some threshold Dc — we see numerically that the 
transmission of the wave decays exponentially with D 
(see Fig. |3]). As a result, the wave propagation depends 
on the number of defects in the chain in this regime. 
This exponential decay arises in the time-evolution of 
the wave amplitude (as the wave propagates through one 
defect after another). The decay of the wave, in turn, 
induces background excitations. Due to the presence of 
the defects, we observe states that seem to be somewhat 
reminiscent of the phenomenon of Anderson localization. 
For a detailed discussion of acoustic analogs of Anderson 
localization and excitation of the pertinent modes in lin- 
ear disordered acoustic chains, see Ref. [32|. We discuss 
the similarities and differences between the phenomena 
that we observe and those known from studies of Ander- 
son localization in more detail in Section I VIII Our ob- 
servation that disorder gradually alters the properties of 
wave propagation is to be expected, given that this also 
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FIG. 2: [Color online] (a) Averaged experimental and numer- 
ical plots of force versus time (as measured by a sensor in 
the 20tli cell of a steel:Al chain) for configurations with a low 
{D = 0.1) and a high {D — 0.7) level of disorder. We normal- 
ize the force based on the peak value for the perfectly ordered 
system (for which D = 0). (b) Normalized peak transmit- 
ted force (sensor in 20th cell; steebAl chain) as a function of 
disorder. Each color/shape represents a different configura- 
tion of equivalent disorder, and each marker corresponds to 
a separate experimental run. The solid curve gives the me- 
dian transmitted force value for each level of disorder, and 
the squares (showing the same qualitative behavior) give the 
results of the numerical runs that correspond to the median 
experimental configurations. 



happens in linear settings. However, for high disorder 
{D > Dc), the response of the chain becomes indepen- 
dent of the level of heterogeneity of the system. This is a 
fascinating property of highly nonlinear chains; beyond 
a critical level of disorder Dc, the wave behavior is in- 
sensitive to the level of heterogeneity in the chain. As a 
result, we see (based on our numerical observations) that 
the peak force of the transmitted wave is well described 
by 

Ft=Fo^ {D->Dc), (3) 

where Ft is the peak force of the transmitted wave, the 
original wave's peak force is -Fb = 1 by normaUzation 
with respect to the maximum force in the perfectly peri- 
odic chain (which has D ~ Q), /i = 3/5 is universal, and 
a « 28 and /3 « 4.4 arc constants whose values depend 
on the particle geometries (i.e., their shape) and mate- 
rial properties in the chain. We measured the values of 
a and /? using numerical fitting for our configuration — a 
large steehsmall steel diatomic chain (the mass ratio is 
mi/m2 = 0.25). We show the force transition in both 
regimes for a large steehsmall steel chain with N = 100 
cells in Fig. [3] (In illustrating our results, we focus 
the discussion in the remainder of this paper on large 
steehsmall steel chains. Our results do not depend on 
this choice, except for particular numerical values when 
indicated.) In the inset of the figure, we show the depen- 
dence of the transmitted force amplitude Ft on the chain 
length N when D ':$> Dc- Namely, Ft decays as a power 
law in N with a coefRcient n that reflects the nonlinear 
relationship between the force and displacement and the 
ID nature of the configuration. The value of the power 
arises from the energy-force scaling oi E p5/3 
the Hertzian contact law (i.e., F ^ S^^'^) and, for exam- 
ple, we would instead obtain fi = 1/2 for a linear contact 
law (for which E ^ F"^). In the next subsection, we will 
justify the effect of the chain length when D <^ Dc- In 
Section lVl CI we discuss the high-disorder regime and the 
value of ^. 

The presence of two regimes with very different dy- 
namics in the aggregate response of the disordered chain 
raises the question of the origin of the transition between 
them. Below we analyze the spatio-temporal structures 
of the propagating waves in both regimes and reveal their 
physical origins. 

B. Low-Disorder Regime: Solitary Waves with 
Decaying Energy 

For concreteness, we frame our discussion using a 
steehsteel chain with 7711/7712 = 0.25 and TV = 100 cells 
(as in Fig. [3]), so that Dc ~ 0.35. When D < D^ the 
propagating wave consists of a leading pulse that resem- 
bles a solitary wave [see Fig. HKa)]. Indeed, the width 
and the shape of the leading pulse are similar to those 
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FIG. 3: [Color online] Semi-log plot of the transmitted am- 
plitude as a function of D for a long spin chain (A*' = 100 
steehsteel dimer cells, with particle sizes chosen to give mass 
ratio mi/m2 = 0.25). Straight lines represent fits of the data 
to Eq. Q. In the inset, we show in logarithmic coordinates 
the effect of the chain length A'' (ranging from = 25 to 
A'' = 300 dimer cells) on the transmitted force amplitude in 
the high-disorder regime D ^ D^. We observe a power-law 
dependence with exponent 0.62 obtained from a fit of the 
data (in agreement with the theoretically predicted value of 
/. = 3/5). 



for (fully-ordered) diatomic chains with D = [10| . How- 
ever, because of the defects, the amplitude of the solitary 
wave decays as it propagates through the chain. The in- 
set shows the peak force Ff of the pulse as a function 
of its position Xf/d in the chain. This force is averaged 
over various configurations of the chain with the same 
level of disorder (i.e., with the same value of D). Note 
that Xf denotes the position of the peak force [so that 
Ff = F{xf)\, and d = R -\- r (where R is the radius of 
the first type of particle and r is the radius of the sec- 
ond type of particle) denotes the length of the cell. As 
the solitary wave propagates, its peak amplitude decays 
exponentially: 



F{xf) = FQe 



(4) 



The length scale ^, which we measure from a fit to 
the numerical results, is approximately 92d for the large 
steehsmall steel configuration in Fig. 21 The value of ^ 
depends on D and on the geometries and material prop- 
erties of the particles. We show the exponential decay 
in the inset of Fig. ID^a). One might be inclined to think 
of the parameter ^ as an analog of Anderson localization 
length in disordered linear settings. A key difference with 
Anderson localization, however, is that the exponential 
localization that we observe occurs in time (as the wave 
propagates through one defect after another), whereas 
localization of Anderson modes occurs in space. Keep- 



ing this key difference in mind, we nevertheless discuss 
a possible analogy with linear systems in more detail in 
Section IVHI 

The decrease of the wave amplitude is a consequence 
of the scattering of the wave as it crosses a defect (i.e., a 
flipped spin) in the chain. To illustrate this process, we 
show the interaction of the solitary wave with an isolated 
defect in Fig. IDJb) using a spatio-temporal depiction of 
the wave as it propagates through a chain with a single 
defect (so that D = 0.02). Observe the two additional 
waves (in addition to the primary pulse) that are emitted 
at the location of the collision between the incident wave 
and the defect. 

The generation of additional waves in the chain results 
in a transfer of energy from the leading pulse to the tail 
[soj . As proposed in Ref. [3l| in the context of weakly 
nonlinear waves in disordered media, this process can 
be described by introducing a transmission coefficient T 
across a single defect. (The coefficient depends on the 
geometries and material properties of the particles.) This 
allows us to express the energy of the leading pulse after 
propagation through a defect as 



E. 



after 



TEi 



before ■ 



(5) 



To measure the microscopic parameter T, we use di- 
atomic chains with TV = 100 cells that contain a single 
defect, and wc measure the amplitude of the wave's force 
profile F both right before the defect and right after it. 
For our calculation, we exploit the fact that the energy 
in the chain scales with force as E p^/^ and that the 
force scales with the displacement S as F ^ S^^^ for beads 
with Hertzian contact [1[ . The transmission coefficient is 

then given by T = [J^after /^before! 

We illustrate this diagnostic in Fig. [5] for a steel: steel 
chain with mass ratio mi/m2 = 0.1. Importantly, the 
transmission coefficient does not depend on the energy 
of the incident wave, in contrast to what has been ob- 



served in investigations of weakly nonlinear waves [31 
For example, we find that T « 0.882 for the steehsteel 
chain of Fig. [5] (which has a mass ratio of mi/m2 = 0.25) 
and that T « 0.84 for the stcehAl chain. To illustrate the 
effect of dimer composition on the value of T, we plot T 
as a function of the mass ratio of the beads in a steel: steel 
chain. We indicate these values with red circles in Fig. [HI 
Note that when raxlmi ~ 1, it is necessarily also true 
that r = 1, as this corresponds to the propagation of a 
solitary wave in a fully ordered (F) = 0) diatomic chain. 

The aforementioned understanding of the effect of a 
single defect allows us to consider the full granular chain 
in the low-disorder regime. Assuming that each defect 
interacts with the wave independently, we can derive an 
expression for the amplitude of the transmitted force 
that results from the successive scattering of the inci- 
dent wave from each individual defect. For a chain with 
disorder D, the expected number of defects that are en- 
countered by the wave after propagating a distance Xf 
is given by 2^0. (Note that obtaining this formula re- 
quires averaging over multiple chain configurations with 
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FIG. 4: [Color online] (a) Wave structure in the low-disorder 
regime (D < Dc) for D — 0.28 and fixed time. We observe a 
clean leading pulse with a shape that is similar to the solitary 
wave that propagates in perfectly ordered chains (which have 
D = 0). In the inset, we show the exponential decrease [see 
Eq. (U)] of the amplitude of the force as the wave propagates 
through the chain. This decrease is induced by the scattering 
of the solitary wave by the defects in the chain. As discussed 
in the text, x/d gives the coordinate normalized relative to 
the dimer size (d is equal to the sum of the radii of the two 
particles), (b) Spatio-temporal depiction of the scattering of 
the wave due to its propagation past an isolated defect (for a 
chain with a single defect, so that D = 0.02). The quantity 
To « 1.7 ms denotes the time that it takes for the wave to 
travel through a perfectly ordered chain {D — 0). 
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FIG. 5: [Color online] Peak force of the solitary wave as it 
propagates past an isolated defect. Far from the defect, the 
wave has the same properties as a solitary wave in the corre- 
sponding ordered chain. The presence of the defect only re- 
sults in a (local) energy loss for the wave as it passes through 
the defect. We quantify this energy loss using the transmis- 



sion coefficient T — 
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FIG. 6: [Color online] Transmission coefHcient as a function 
of mass ratio for the solitary wave that we observe in the low- 
disorder regime [using Eq. ((5])]. The values that we obtained 
from considering the energy loss when passing through an iso- 
lated defect (see Fig.[5l) show excellent quantitative agreement 
with those that we obtained from the aggregate transmission 
properties of the elastic spin chain [using Eqs. ((3l[7ll], which 
we derived from our model of wave propagation in the low- 
disorder regime. 
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the same disorder value.) The pulse energy is reduced to 
E ^ E{t = 0)T^°. The relation E pi implies that 



F(a;/) =Foe- 



lOd 



(6) 



where -Fb = 1 by normalization. This agrees with the 
wave structure in Eq. ([4]) that we obtained numerically. 
Using T « 0.882 and D ?» 0.28, which was the case for 
the example in Fig. UJ we obtain ^ « 95d from in 
excellent agreement with the fit of the numerical data in 
the inset in Fig. lUJa). 

We now consider wave propagation along the entire 
chain by taking Xf = Nd. This yields the exponential 
decay of Eq. (|3]), with parameter 
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3iog(i/r)- 



(7) 



Equation (|3]) can thereby be used to describe the de- 
crease of the transmitted amplitude for D <^ Dc- The 
value a w 27 for the steel: steel chain with mass ratio 
mi/m2 = 0.25 that we obtained from the transmission 
coefficient T « 0.882 using this expression matches the 
value that we measured directly in numerical simulations. 
To further test our theoretical results, we compare in 
Fig. ini the value of T predicted by the exponential de- 
crease of Ft{D) using Eq. ([3|) with the value that we 
measured from the interaction of the wave with a sin- 
gle defect. These results agree quantitatively within 1%, 
confirming that the defects can be treated separately 
when considering the transmitted force, provided that 



C. High-Disorder Regime: Delocalization of the 
Waves 



We saw earlier that if the level of disorder in the chain 
is low (e.g., if D <C Dc), then the disorder D governs 
the wave transmission through the chain by fixing the 
number of defects that scatter its energy during propa- 
gation. However, when there are too many defects, the 
amplitude of the transmitted wave becomes robust to the 
configuration of the clastic spin chain and no longer de- 
pends on D (see Fig. [5]). A detailed analysis of the wave 
structure in this high-disorder regime reveals that the 
spatio-temporal evolution of the wave is now indepen- 
dent of the level of disorder, as the system has reached a 
level of maximum scattering by the chain heterogeneities 
and additional defects do not further delocalize the prop- 
agating wave. 

In Fig. m we show the basic process of wave delocal- 
ization in highly disordered chains. To obtain each curve 
in this figure, we fix the disorder level D > Dc and ran- 
domly assign the positions of the N x D/2 defects in that 
chain. We then measure the force amplitude as the wave 




FIG. 7: [Color online] Spatio-temporal structure of the wave 
in the high-disorder regime {D — 0.70). As in Fig. |4l the 
quantity to 1.7 ms denotes the time that it takes for the 
wave to travel through the perfectly ordered chain. 



propagates through the chain and plot the mean 
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where F\x) gives the force amplitude of the jih con- 
figuration. As the wave propagates, the force amplitude 
Ff = F{xf) of its front decays and the wave broadens 
in space so that its energy spreads throughout all (con- 
tinuous) positions x G [0, a;/]. To quantitatively describe 
the spatial and temporal dynamics of the wave struc- 
ture, we average over the different configurations of the 
elastic spin chain with disorder parameter D in order to 
determine the primary wave dynamics at a fixed level of 
heterogeneity. Because the aggregate transmission of the 
wave is robust to the level of heterogeneity in the chain, 
we find that the averaged spatio-temporal wave structure 
is also independent of D. It is given by 



F{x,t)^H{xf~x)Ffe 



(8) 



where H is the Hcaviside function, xj is the position of 
the front of the delocalized wave, and Ff = F{xf) is 
its corresponding amplitude. As illustrated in Fig. E 
the decay of the force amplitude of the wave front is 
well described by a power law Ff ~ xJ°-^\ in excellent 
agreement with the theoretical prediction of /i = 3/5, as 
it propagates through the chain. In contrast to the wave 
dynamics in the low-disorder regime, the wave scattering 
in the high-disorder regime is independent of the disorder 
parameter D. 

We can capture the main properties of the wave in the 
high-disorder regime by exploiting the fact that the scat- 
tering has reached its maximum level in the chain and 
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Fig. [H From a mathematical perspective, this imphes 
that we can write the wave structure as 
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FIG. 8: [Color online] Force amplitude of the leading front 
of the delocalized wave as a function of its position Xf in the 
chain. For each curve, we average over 10 configurations of 
the chain with disorder level D > Dc- As the wave propagates 
and after a short transient, the force amplitude decays as a 
power law with a coefficient of about 0.63, which is close to 
the theoretical prediction of = 3/5. 



approximating the resulting wave structure by a highly 
delocalized wave in which the elastic energy is spread out 
evenly among the dimer cells for which x G [0,a;/]. (The 
only cells that are excluded are those that the wave has 
not yet had enough time to reach.) In this situation, the 
scattering in the chain has reached its maximum level 
because the wave is now fully delocalized. Additional 
defects do not cause further delocalization, so the wave 
properties are robust in this regime to the specific value 
of D. Such a stationary spatial structure of the wave 

yields an energy per cell of -^yj^ when the front reaches 

Xf. Using the relation F ~ E'^/^ between the peak force 



and the wave energy yields F 



-3/5 



'see the discus- 



sion below), which agrees with the numerical observa- 
tions of Fig. [5] If we consider such stationarity over the 
full chain — i.e., when the front reaches the end of the 
chain, so that Xf = Nd — we obtain the force amplitude 
of the transmitted wave through the chain. This gives 
a plateau whose level is F ~ A^~3/5 jggj, jgjjj^ g^. 
cellent agreement with the direct numerical simulations 
shown in the inset in Fig. [3l 

Although one can observe behavior like that shown in 
Fig. [5] when there is energy equipartition, such a mecha- 
nism is not necessary. Indeed, as one can see from Eq. ([8|), 
we do not have equipartition. Instead, we observe that 
the wave structure F(x,t) remains stationary while ex- 
panding, which is necessary to obtain the features in 



F{x,t)=Ff{xfy 



Xf 



(9) 



where u is the spatial structure of the wave, which ex- 
pands as the position of the front xj increases. The am- 
plitude of the wave and its evolution during the propa- 
gation is given by Ff{xf). By energy conservation, the 
integral of the energy density F{x, t)^/'^ over the portion 
of chain where the wave is non-zero — i.e., over the inter- 
val [0, Xf] — must remain constant. It is thus independent 
of X/ , which implies that the amplitude of the wave de- 
cays as Ff(xf) Xj^^^ . In contrast, in the solitary-wave 
regime, the structure of the wave is not stationary, and 
energy is transferred from the leading pulse to the tail of 
the wave, so that the wave structure cannot be written 
in the form In the high-disorder regime, the fluxes of 
energy from the leading part of the wave to the tail and 
those from the tail to the leading part of the wave bal- 
ance each other, and the wave structure docs not change 
aside from the dilation in Eq. ([9]). 



D. Finite-Size Effects 

From our understanding of the low-disorder and high- 
disorder regimes, we can now analyze the effect of the 
chain length on the transition from solitary wave propa- 
gation (for D ^ Dc) to highly delocalized waves {D 3> 
Dc). We have used the force amplitude Ft of the wave 
as a function of the level of disorder D in the chain to 
characterize this transition: The force Ft first decreases 
exponentially with D and then saturates to a plateau for 
D ^ Dc- This provides a natural way to define the crit- 
ical level of disorder Dc as the intersection between both 
regimes [see Fig. (Sla)]- Equating the two expressions in 
Eqs. which give the force of the incident wave as a 
function of the chain length in the two regimes, we obtain 
the critical defect density 



Dc 



N 



[fLiogiN)~iogm 



(10) 



where a is expressed in Eq. ([7]) as a function of the trans- 
mission coefficient of the wave through one defect and 
/i = 3/5 is a universal exponent. As the chain length 
increases, the transition occurs for smaller values of Dc, 
and the threshold tends to as the chain becomes in- 
finitely long. This result can be interpreted as follows: It 
is the number of defects that the wave encounters rather 
than the defect density that determines the type of wave 
propagation (i.e., whether one observes a solitary pulse 
or a delocalized one). As a result, for a very long chain, 
even a low defect density is able to generate delocaliza- 
tion. When the force amplitude of the leading pulse is 
comparable to that of a rcfiected wave, we are in the 
delocalizcd-wave regime. This occurs when the ampli- 
tude of the leading pulse is below a critical level, which 
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happens once it has encountered a given number of de- 
fects (recaUing that the amphtude decays by a factor of 
1 — r each time the wave encounters a defect). 



VII. DISCUSSION 

Let us now compare and contrast our observations of 
the wave properties that we have seen in highly non- 
hnear chains with those reported previously for linear 
and weakly nonlinear settings (see Ref. [32] for a re- 
view). In particular, it is interesting to consider the wave- 
scattering mechanisms that we have observed in elastic 
spin chains versus the process of Anderson localization 
that can occur in linear and weakly nonlinear systems 

MM- 

We have observed that the heterogeneities present in 
the granular chain alter the wave propagation, so that 
in the limit of a very long chain, even a low defect den- 
sity is sufficient to prevent wave propagation through the 
whole system. This intuitive effect is also observed in lin- 
ear disordered chains through the localization in space of 
the wave close to the excited part of the chain for ar- 
bitrary levels of disorder, as ensured by Furstenbcrg's 
theorem [s^. (Here, we use the terms "localization" and 
"localized wave" to contrast these observations with the 
extended state of the waves that arise in periodic and 
perfectly ordered linear chains.) A signature of Ander- 
son localization in linear disordered settings is the ex- 
ponential decay of the wave amplitude as it propagates 
through the chain [s^ [s^. This defines the Anderson 
localization length involved in such an evolution law. As 
we have reported in this paper, elastic spin chains ex- 
hibit exponential amplitude decay [see Eq. (|4])] for low 
levels of disorder {D <^ Dc). It is then perhaps tempting 
to consider the localization length ^ defined in Eq. ([6]) 
as a nonlinear counterpart of the Anderson localization 
length for highly nonlinear settings. However, the reader 
is cautioned towards such a connection, as the decay of 
the amplitude is connected to the wave propagation over 
time (as additional defects arc encountered) and does not 
really relate to the structure of the spatial profile of the 
wave. Furthermore, recall (see , for example, the discus- 
sion of Section 4.2 in Ref. [33) that FPU-type systems 
do not admit a genuine form of Anderson localization (at 
least for small frequencies) unless a harmonic on-site po- 
tential is added to each node in the lattice. As another 
point of comparison, we have provided in this paper an 
expression for ^ as a function of the defect density D 
and the scattering characteristics of one isolated defect 
T, which is also possible for linear disordered chains in 
some specific situations [13, . 

Localization in highly nonlinear systems can be funda- 
mentally different than localization in linear and weakly 
nonlinear systems. In particular, we observe in granular 
chains a fascinating property (which is a central result 
of our study): Beyond a critical level of disorder, the lo- 
calization length saturates, so that the wave dynamics 



remains insensitive to the level of heterogeneity in the 
chain in this regime. In addition, the decay of the wave 
amplitude does not follow an exponential relation but 
instead follows a power-law relation (see Fig. [5|), with a 
power fixed by the type of nonlinear interaction between 
nearest neighbors in the chain. This remarkable and un- 
expected effect arises from the balance between the flux 
of energy from the leading pulse to the tail of the wave 
and the flux in the opposite direction. It results in a 
stationary wave structure that, combined with an argu- 
ment based on energy conservation in the chain (see Sec- 
tion (VTC]), suggests the possibility of power- law decay 
(in time) of the wave amplitude for sufficiently high lev- 
els of disorder. How and when this condition of flux bal- 
ance is achieved still remains to be established. However, 
our investigation of the wave behavior in both the low- 
and high-disorder regimes gives an expression of the crit- 
ical level Dc of disorder corresponding to the transition 
from exponential to power-law decay of wave amplitudes 
[see Eq. ([TU| ]. It is noteworthy that for sufficiently long 
chains, the system experiences the second type of decay 
(because ^ \og{N)/N). Consequently, this regime 
is expected to dominate for long times and large system 
sizes. 



VIII. CONCLUSIONS 

We have introduced the concept of elastic spin chains, 
which wc constructed using randomly-oriented arrange- 
ments of dimer cells in one-dimensional granular crystals. 
We quantified the level D of disorder in the chain by 
counting the number of spins with the same orientation. 
In order to characterize the wave propagation proper- 
ties, we calculated ensemble averages of relevant quan- 
tities over different chain configurations with the same 
value of D — i.e., configurations with the same number 
of up and down spins. This allowed us to quantitatively 
characterize the effect of disorder on the propagation of 
highly nonlinear pulses. For low levels of disorder, we ob- 
served the propagation of a solitary wave with exponen- 
tially decaying amplitude. We observed amplitude decay 
in time (as the wave propagates through one defect after 
another). The "localization length" characterizing the 
exponential amplitude decay in the low-disorder nonlin- 
ear setting can be expressed as a function of the level of 
disorder D and the scattering properties of an isolated 
defect (i.e., when one dimer has a different spin from 
all of the other dimers). As a result, the transmission 
coefficient through the entire chain decays exponentially 
with the level of disorder D. Beyond a critical level of 
disorder Dc, this transmission coefficient becomes inde- 
pendent of D and saturates at a constant value. In this 
regime, in which the wave behavior remains insensitive to 
the level of disorder, the wave's spatial structure is much 
more extended. Moreover, the propagation of the wave 
is no longer governed by the propagation of an isolated 
solitary wave, as it was the low disorder regime, but it 
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instead resembles a train of smaller waves. In addition, 
the decay of the wave amplitude follows a power law in- 
stead of an exponential law, with the power fixed by the 
type of nonlinear interaction between nearest neighbors 
in the chain. This effect appears to arise from a balance 
between the flux of energy going from the leading pulse 
to the tail of the wave and the flux going in the opposite 
direction. How and when this condition of flux balance is 
achieved remains to be established, and additional stud- 
ies of localization phenomena in highly nonlinear settings 
should lead to interesting insights into the dynamics of 
disordered nonlinear systems. 
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